Superconductivity in electron-doped cuprates: 
Gap shape change and symmetry crossover with doping 
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The Kohn-Luttinger mechanism for superconductivity is investigated in a model for the electron 
doped cuprates. The symmetry of the order parameter of the superconducting phase is determined 
as a function of the geometry of the Fermi surface together with the structure of the electron-hole 
susceptibility. It is found to remain d x 2_ y i wave within a large doping range. The shape of the gap 
anisotropy evolves with doping, with the maximum gap moving away from (n, 0), in good agreement 
with recent experiments. As the shift of the maximum increases, a crossover to d^-symmetry is 
found. 
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I. INTRODUCTION 

The single band Hubbard model with longer range hop- 
pings - (t-f) Hubbard model - is being widely accepted 
as a minimal model to describe the physics of the high- 
T c cuprates [1]. It gives most of the qualitative aspects 
of the phase diagram: antiferromagnetism near half fill- 
ing, superconductivity, pseudogap, and striped phases. 
The greatest effort has previously been centered in un- 
derstanding hole-doped compounds as they possess the 
highest available T c . The complication of the phase di- 
agram in the underdoped to optimally doped regimes, 
partially due to the inhomogeneous structures and the 
proximity of the Van Hove singularity, has recently re- 
newed interest in electron doped compounds [2]. 

Although the estimated values of the Hubbard inter- 
action U lie in the intermediate regime, most of the pre- 
vious features can be obtained in the weak coupling limit 
which has the advantage of often permitting analytical 
computations where the physics is more transparent. 

Superconductivity has been obtained in the purely re- 
pulsive model with quantum Monte Carlo [3] and other 
numerical methods [4], with the antiferromagnetic fluc- 
tuation exchange approximation [5-7], and with Renor- 
malization Group (RG) methods in the proximity of the 
Van Hove singularity [8]. At present there is no gen- 
eral agreement on the mechanism for superconductivity. 
As for the symmetry of the superconducting order pa- 
rameter, it is generally assumed to be d-wave in the hole 
doped cuprates and remains unclear in the electron doped 
compounds [9]. 

One of the most physically appealing mechanisms for 
obtaining a pairing instability from repulsive interactions 
goes back to Kohn and Luttinger (KL) [10] who demon- 
strated the instability of a three dimensional isotropic 
Fermi system towards pairing due to the spatial modu- 
lation of the effective interaction at a wavevector of 2kp. 
The KL mechanism has been extended to other electron 



systems and dimensions [11-13] and has been extensively 
studied in the 2D Hubbard model [14-17]. The Kohn- 
Luttinger mechanism has very recently been reanalyzed 
in [18]. 

In the study of the instabilities of an electron system 
a principal role is played by the electron-hole suscepti- 
bility w) whose imaginary part measures the den- 
sity of electron-hole pairs with energy to and whose real 
part renormalizes the scattering amplitude. Scattering 
processes which involve opposite points of the Fermi sur- 
face can be enhanced by its special geometry (nesting 
and Van Hove singularities are extreme examples), or by 
other physical features. 

In ref. [19], a scaling analysis was used to study the 
pairing instabilities of a general Fermi surface in the 2D 
square lattice as a function of its geometry. It was shown 
that the curvature of the Fermi line modulates the effec- 
tive interaction in the BCS channel in such a way that 
different harmonics scale as different powers of the scal- 
ing parameter. As the latter goes to zero some harmonics 
become negative giving rise to a KL superconductivity in 
the given channel. 

In this paper we will perform a KL analysis of elec- 
tron doped cuprates along the lines of ref. [19] based on 
special features of the susceptibility. Recent experiments 
on electron doped cuprates [9] propose a change in the 
symmetry of the superconducting order parameter from 
d wave below and around optimal doping, to s wave in 
the overdoped regime. Moreover, even when the gap has 
overall d-wave symmetry, its angle-dependence evolves 
with doping, picking up substantial harmonic content 
for hole underdoping, and there is some evidence [20] 
that the peak shifts away from (tt, 0) in electron-doped 
cuprates. This issue will also studied in the paper. 

The organization is as follows: In sect. II we set the 
model, review the main arguments of the scaling analysis 
of ref. [19] and establish the d-wave nature of the super- 
conducting phase. In sect. Ill we analyze the structure 
of the spin susceptibility to be used in the calculation. 
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Sect. IV is devoted to the evolution of the symmetry of 
the order parameter with doping followed by our conclu- 
sions. 

II. THE MODEL AND KOHN-LUTTINGER 
MECHANISM 

The Fermi surface of the electron compounds for the 
doping values of interest has a general rounded shape 
centered at (jr, n) with flatter portions in the diagonal 
directions. The the t-t'-Hubbard model with negative 
values of t' is the simplest model that reproduces the 
observed feature. The dispersion relation is given by 

e(k) — — 2i(cos k x + cos k v ) — At' cos k x cos k y . (1) 

For definiteness we assume t = 0.326eV, t' jt = —0.276 
in our calculations [21-24]. 

The Kohn-Luttinger mechanism is closely related to 
Friedel oscillations. It is well known that, due to the 
sharp cutoff of the electron distribution in /c-space at the 
Fermi level, impurity potentials in a metal do not fall 
off monotonically but have a superposed (Friedel) oscil- 
lation. Kohn and Luttinger showed that a similar oscilla- 
tion arises in the electron-electron interaction, leading to 
an attractive interaction between two electrons separated 
by the right distance - the position of the first Friedel 
minimum. In turn, the attractive interaction can lead to 
a superconducting instability. 

The calculation can readily be formulated in renormal- 
ization group language [25], based on the fact that the 
effective coupling constants of a given hamiltonian (ver- 
tex functions) acquire an energy-momentum dependence 
upon renormalization behaving like effective potentials. 
In the simplest Fermi liquid model (gas of electrons with 
spherical Fermi surface and short range four fermi inter- 
actions), only the forward and BCS channels get renor- 
malized [25]. The standard KL mechanism occurs when 
the effective BCS vertex at a given momentum i2kp) os- 
cillates in such a way that some of its Fourier components 
becomes negative. The system then undergoes a super- 
conducting transition. The symmetry of the supercon- 
ducting order parameter can be found by expanding the 
potential in eigenfunctions of the symmetry of the model 
(spherical harmonics in the spherical case) and finding 
the lowest negative eigenvalue. 

This is the analysis that we will follow in sect IV of 
the paper adapted to case the Fermi surface given by the 
contour lines of eq. (1). 

The KL mechanism for the rounded Fermi surface 
corresponding to electron doped cuprates was analyzed 
in [19] and shown to induce a pairing instability with 
d x 2_ y 2 symmetry without special features of the suscep- 
tibility. Based on very simple scaling arguments, it was 
shown that the electron susceptibility is proportional to 



(l// 2 ) where / is the curvature of the Fermi surface and 
gets modulated by it. Hence it has maxima for the scat- 
tering vectors joining two opposite points of the Fermi 
surface in the (tt/2,tt/2) direction and equivalent points 
where the curvature reaches its minimum value. It has 
minima in the zero and equivalent directions. This situ- 
ation corresponds to d x i_ y i symmetry This analysis is 
independent on the specific form of the spin susceptibil- 
ity. The detailed study of the susceptibility in the next 
section reinforces the symmetry arguments of [19] and 
the d-wave character of the instability. 

III. THE STRUCTURE OF THE SPIN 
SUSCEPTIBILITY 

Recent angle-resolved photocmission spectroscopy 
(ARPES) of electron-doped Nd 2 _ 2 Ce ;E Cu04± ( 5 (NCCO) 
[2] found a smooth evolution of the band dispersion with 
doping, from half filling to optimal doping, which could 
be interpreted [26] in terms of the gradual filling of the 
upper Hubbard band. Doping gradually reduces the 
Mott gap that closes near optimal doping. This dop- 
ing also falls close to another interesting point [27] where 
the displaced Fermi surface is tangent to the original one 
and the susceptibility x at Q — (n, 7r) starts to drop pre- 
cipitously. 

Near this point, the magnetic susceptibility has the 
form of a nearly flat-topped plateau in momentum space, 
with sharp falloff away from the plateau. The plateau 
is defined by the presence of points at which the Fermi 
surface (FS) overlaps the FS image shifted by Q. If the 
FS image is shifted away from Q by an additional q ', 
then for some critical value q 1 = q c , one or more of the 
overlaps ceases to exist as the two FS's pull apart. The 
value of q c , which defines the plateau boundary in a given 
direction, satisfies 

-2i(sin( 9x /2) + sm(q y /2)) - At' sm(q x /2) sm{q y /2) = M , 

(2) 

where \i is the chemical potential. 

The plateau exists in the doping range > /i > 
I^vhs = 4i', where hvhs is the doping of the Van Hove 
singularity. When /j, = 0, the width of the plateau shrinks 
to zero. The plateau approximately satisfies the form 
X = A — BQ(q' — q c )\/q' — q c , with a step function. 
As the width of the plateau vanishes q c — ► 0, the two 
square root terms merge at a single point. The result- 
ing susceptibility is depicted in Fig.l. We find that the 
angle dependence of the resulting superconducting gap 
function changes dramatically depending on whether the 
doping lies on or off of this susceptibility plateau (leading 
ultimately to changes in the gap symmetry). 
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FIG. 1. Bare susceptibility at the termination of the 
plateau described in the text, [i — 0. Brillouin zone points 
are T = (0,0), X = (tt,0), S = (tv,ty). 

At T=0, the susceptibility \q can be written 

d 2 k 1 



FIG. 2. (a) Calculated susceptibility xil) f° r several values 
of overlap ks (from right to left, kg/q = 0.5, 0.25, 0.05, and 0, 
the last corresponding to the termination of hot spot overlap), 
(b) Schematic of Fermi surfaces, defining kg, k± and kit. 

The FS centered at (tt, tt) and the Q-shifted FS are 
illustrated schematically in Fig. 2b. The region of inte- 
gration is over the part of the upper FS in Fig. 2b not 
overlapped by the lower (Q-shifted) FS, and k± ranges 
from zero at the apex of the wedge to k c = kp — kg at 
the middle of the upper FS, where Hf is the radius of the 
FS and kg is the overlap parameter defined in Fig. 2b. 
Assuming kg << kp and keeping only the lowest order 
contributions, 



d 2 k 



1 



dk 



dk 



(2tt) 2 



2kp + ykpkg In 



y/kg - kj_ 
k± 1 
1-/3, 



1 + /3 1 



(6) 



2tJ (2 tt) 2 (cos (k x + q x ) + cos (k y + q y ) — cos k x — cos fc y )with^/3 = y/kg /kp. For the g-dependent susceptibil- 
ity, let q = Q + q'. Then the FS is shifted by q' , or 
kg — ► kg + q' /2 for the surface shown, kg — > kg — q'/2 
for the FS at (-tt, -tt), so \ q « h s +q'/2 + h e - q '/2- 
There is one correction. For q' > 2kg, the two FSs 
no longer overlap, and the integral runs over the full 
half-circle, but k± is measured from half way between 



The Fermi functions limit the integral to a sum of ap- 
proximately wedge-shaped areas. We begin at the an- 
tiferromagnetic wave- vector, q = Q = (tt,tt). Letting 
ki = n/2 + k[, i = x,y, then to lowest order the energy 
becomes 



2v / 2tfc'i 



2t'k\f 



(4) 



with fen and k± the momenta parallel and perpendicular 
to the zone diagonal (magnetic Brillouin zone boundary) . 
Linearizing the energy denominator, Ae oc kj_, indepen- 
dent Of fell 



the two FSs 

7 = VW 2 



l k s -q'/2 



dkidk 



k± 



4ttH 



(5) 




0.5 1 1.5 q 2 




2k F [l - 7 tan" 1 1/7], with 

2kg) /2kF- Figure 2a shows the resulting 
susceptibilities for several values of kg. The calculation 
explains the very flat top with weak positive curvature 
and the sharp falloff at q 1 = 2kg. For increased electron 
doping, the hole-like Fermi surface shrinks, kg — ► 0, and 
the plateau width shrinks to zero. Near the plateau edge 
X ~ 1 — 7T7/2 varies as ^fq 1 . The cusp-like susceptibil- 
ity in Fig. 2a at the point where the plateau terminates 
corresponds to the xo peak in Fig. 1 at S = (tt,tt). 
With the assumed parameter values, this happens when 
fj, = 0, corresponding to an electron doping x ~ —0.19; 
on the hole doped side, the plateau terminates when 
H = At' = -0.359ey, at a hole doping, x = 0.24. 



IV. COUPLING CONSTANT CALCULATIONS 

Calculation of the magnitude of the superconducting 
transition gap is beyond the scope of the present paper: 
this requires a better understanding of (a) the proper 
choice of susceptibility, (b) incorporation of the frequency 
dependence of the susceptibility, (c) proper accounting of 
the competition with magnetic ordering, and (d) solution 
of the resulting (generalized) Eliashberg equations. 

The most negative coupling constant determines the 
dominant gap symmetry, while the corresponding eigen- 
function gives the angle dependence of the gap. 
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In this section we investigate the possible change in 
the symmetry of the superconducting order parameter 
and how its angle-dependence evolves with doping [20]. 

The pairing coupling constant in a given symmetry 
sector is given by the matrix [6, 17] 



(2 



i rdk r ^^ (k)k / )Aan(k)Aam(k0 
J Vk J v k t 



(J) 

where V(k,k') = U + f/ 2 x(k + k'), and A an (k) is the 
normalized ( f — A Qm A Q „(k) = S mn ) weight function ex- 
panded in terms of the irreducible representations h an of 
the symmetry group. We approximate the CuC>2 plane by 
a square lattice, in which case the appropriate symmetry 
group is D4, for which there are four singlet and one dou- 
blet representations. These representations define gap 
symmetry sectors: while Eq. 7 can mix basis functions 
within a given sector, it does not mix functions between 
sectors. The four singlet sectors are labelled according 
to their lowest basis functions, as s, d (for d x 2- y z), d xy , 
and g, while the doublet sector is labelled p - actually, 
there are two subsectors, p x and p y which do not mix, 
but have degenerate eigenvalues. The corresponding ba- 
sis functions are (with n ranging from to 00) 



h s ,n(<t>) 


= cos[4n0] 






hd,n(4>) 


= cos[(4n - 


-2)0] 






= sin[(4n H 


-2)0] 




hg, n {4>) 


= sin[4(n H 


-1)0] 






= sin[(2ro H 


-m 






= cos[(2n -I 


-i)0]. 


(8) 



The above functions describe the orbital symmetry - the 
singlet (doublet) representations corresponding to spin 
singlets (triplets). Below, we evaluate the lowest eigen- 
value for each symmetry sector, and label the correspond- 
ing eigenfunction by the sector - even though, e.g., the 
lowest s-wave solution is really an 'extended-s' solution, 
with all the /i Sj „, n > 0, orthogonal to the pure s-wave 
solution. 

For fi < /j-vhs, the Fermi surface is an electron-like 
Fermi surface closed about the T point (0,0), and the 
angle must be measured about this point. For fi > 
Hvhs, the topology of the Fermi surface changes to hole- 
like, centered on (tt, tt), and must now be measured 
from this point. 




FIG. 3. Development of coupling constant with increas- 
ing matrix size, N, for p, = — O.leV. Circles = d x 2_ y 2\ 
up-pointing triangles = g; diamonds = p x ; squares = s; 
down-pointing triangles = d xy . 

The maximal superconducting coupling is given by the 
minimal (i.e., maximally negative) eigenvalue of the A 
matrix, Eq. 7. To solve this equation, the matrix was 
cut off at a finite size N x N, with N — 15. This large 
N value was employed to assure adequate convergence 
in all sectors. This is illustrated in Fig. 3. If the N=l 
eigenvalue is already negative, the only change with N 
is a small increase in magnitude (due to level repulsion) . 
But if the N=l eigenvalue is positive, N plays a larger 
role. There are two effects: first, a diagonal matrix ele- 
ment might itself be negative and second, level repulsion 
always pushes the largest and smallest eigenvalues away 
from the mean. For example, for the s-wave sector, all 
the diagonal elements are found to be positive, but level 
repulsion generally leads to a small negative eigenvalue. 
However, for the parameter range studied, this is always 
too weak to be of interest. 
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FIG. 4. (a) Evolution of the superconducting effective cou- 
pling as a function of doping, for a variety of gap symmetries: 
solid line = d x 2_ y 2\ dot-dashed line = g; dotted line = p x ; 



long dashed 
doping. 



s; short dashed = d xy . (b) Blow-up of electron 
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The results of the analysis are shown in fig. 4 for a 
variety of dopings (for electron doping x is considered 
negative), assuming U = 6t. Reasonable A values are 
found for both electron and hole doping, but do not in- 
clude the suppression of superconductivity near half fill- 
ing caused by the magnetic order. While A decreases with 
electron doping, the preferred symmetry remains d x 2_ y 2 
over the doping range of interest. The present calcula- 
tion thus provides no indication for adtos crossover of 
the symmetry. However, at a higher doping, x ~ —0.39, 
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there is a crossover from d x 2_ y 2 to d^j, symmetry; near 
such a crossover, there is likely to be a range of (gap- 
less) d a; 2_ y 2+id 2 ;j, symmetry, which may simulate an s- 
wave gap. Alternatively, the symmetry change may be 
associated with an additional pairing contribution due to 
electron-phonon coupling. It is interesting to see that d- 
wave symmetry is dominant at both sides of half filling, 
as found experimentally. 




FIG. 5. Normalized coupling constants, X/Z, for same data 
as in Fig. 4. 

The Van Hove instability is clearly noticeable as the 
(divergent) peak in the hole-like sector. However, for su- 
perconductivity, the relevant parameter is likely to be 
X/Z, Fig. 5, where Z = 1 + A s and X s is the (0,0) ma- 
trix element in the s-wave sector of Eq. 7 (note that this 
is the only term which includes the linear-in-£7 contri- 
bution to V). This renormalization eliminates the Van 
Hove peak, shifting the largest d-wave gap to a lower hole 
doping. This is quite suggestive of the experimental situ- 
ation, where there appears to be a quantum critical point 
somewhat above optimal hole doping [28] , which may be 
associated with the Van Hove singularity (VHS) [27]. 

We find a striking evolution of the angular dependence 
of the d-wave gap, as shown in Fig. 6. The Kohn- 
Luttinger mechanism leads to significant harmonic ad- 
mixture, which changes as a function of doping. In- 
cluding all harmonics, the gap functions contain exces- 
sive structure, so the following approximation was intro- 
duced, to provide smoother and more robust gap func- 
tions. For each added harmonic order (going from an 
N x N to an (N + 1) x (N+l) matrix), the change in the 
smallest (largest negative) eigenvalue was monitored, and 
if the fractional change was less than some small refer- 
ence value a (typically, a = 0.02 was used), the coupling 
to this harmonic was neglected. This reduced the matrix 
problem from 15x 15 to N x N, where N was generally 
2-4, except in the immediate vicinity of the Van Hove 
singularity, where more harmonics were needed. An ex- 
ample is shown in Fig. 6c, while Figs. 6a, b show only 
smoothed data, except at the VHS, where all harmonics 
are shown in Fig. 6a, and only the dominant N = 7 in 
Fig. 6b. 

For all hole dopings the peak stays close to (tt, 0), 
but harmonic content tends to sharpen the peak and 
flatten A near the nodes, deviating from the simplest 
cos k x — cos k y form. Such flattening has been found in 



both ARPES [29,31] and STM [32] experiments. How- 
ever, the experimental trend is that the harmonic content 
is enhanced in underdoped samples [30] , whereas the cal- 
culated trend is for larger harmonic content to develop 
near the VHS - i.e., with increasing hole doping. [In com- 
paring to experiment, it must be kept in mind that, for 
/i > — 0.3599ey, all angle 4> measurements are centered 
at (tt, tt), with = corresponding to (tt, 0).] 




10 20 30 d) 40 




FIG. 6. (a) Angular dependence of the superconducting 
gap for the d-wave symmetry solution, for several dopings: 
the chemical potentials [hole densities] are n [x] — -0.3599 
[0.247] (thin solid line), -0.35 [0.22] (long dashed line), [-0.19] 
(short dash-dotted line), 0.10 [-0.28] (dotted line), 0.20 [-0.35] 
(dash-dot-dotted line), and 0.3eV [-0.41] (long dash-dotted 
line) . Arrow = datum of Ref. [20] . (b) Continuation to higher 
hole doping, with a [x] = -0.3599 [0.247] (thin solid line), -0.38 
[0.30] (long dashed line), -0.4 [0.33] (short dashed line), -0.45 
[0.41] (dot-dashed line), -0.5 [0.49] (dotted line), and -0.55 
[0.55] (dash-dot-dotted line), (c) gap function for /i = O.SeV, 
comparing a calculation employing all 15 harmonics (dotted 
line) with one involving only the dominant 4 harmonics (solid 
line) . 

As the doping shifts toward electron-like, there is a 
significant shift of the peak position with cf>, away from 
(tt, 0). This result is consistent with recent observations 
by Blumberg, et al. [20] (One should, however, note the 
debate [33,34].) However, the agreement is not quanti- 
tative: a shift of the peak to ~ 15° is found for an ex- 
perimental doping of x = —0.15, whereas the predicted 
doping is —0.32, Fig. 7. Indeed, there appears to be a 
close correlation between the shift of the d-wave peak 
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and the crossover of the symmetry from d x 2_ y 2 to d xy : 
the crossover occurs when the peak has shifted about 
half way to 45° (asterisk In Fig. 7). Note that a sim- 
ilar shift arises on the hole-overdoped side, although a 
second, larger peak remains at = 0. 
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FIG. 7. Shift of d x 2_ y 2 peak from (it, 0) as a function 
of doping (a) and chemical potential (b). The filled circles 
(squares) correspond to changing t' to — 0.12t (— 0A25t); the 
dashed line in (b) is the position of the hot spots. Arrow in 
(a) = datum of Ref. [20]. Asterisk = point of crossover to 
d^y-symmetry. 



Although the calculations are based in a weak coupling 
analysis we believe that they are justified in the electron- 
doped case and that more refined computations will not 
change the general features. 

Note added: Recently, we became aware of a simi- 
lar calculation [35]. Here, a doping-independent nearly- 
antiferromagnetic Fermi liquid susceptibility was intro- 
duced in place of the lowest-order (Kohn-Luttinger) form 
we assumed. A very similar crossover of gap symme- 
try with electron doping was found, with two differences. 
First, the crossover was found to be at the doping at 
which the hot-spot plateau terminated. This is probably 
not very significant, however; the parameters are so dif- 
ferent that this actually corresponds to a higher doping 
(x ~ —0.59) than we found, and much higher than in the 
experiments. More significant is that we find a crossover 
from d x 2_ y 2 to d xy symmetry, while their crossover is to 
p-wave symmetry. We have repeated our calculations us- 
ing the full self-consistent susceptibility [27] , and find for 
this nearly divergent susceptibility that the crossover is 
to a state of either p- or g-wave symmetry. These results 
will be reported in a future publication. 
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This shift of the peak from (it, 0) does not follow the 
position of the hot spots on the Fermi surface (dashed line 
in Fig. 7b), but depends on whether or not the system 
is on the susceptibility plateau. The peak stays close 
to (7T, 0) as long as the chemical potential is on the hot 
spot plateau (/i < 0), then rapidly shifts toward 45°. 
Similarly, for n < — 0.3599eV (off of the other side of the 
plateau, beyond the VHS), a second peak appears and 
shifts to 45° by ~ fi — —0.5, Fig. 6b, by which point the 
gap has crossed over to p-wave symmetry. 

For the electron doping case we have studied how 
this peak shift changes when the plateau width (or t) 
is varied. From Eq. 7 the angle dependence is con- 
trolled by the product of two terms, a weighted dos, 
g(k) = A s .i(k)k 2 (0)/vk, and a weighted susceptibility 
V(k) = Jd0'V(k,k')g(k'). The latter V(k) peaks at 
(-7T, 0) on the plateau, and starts shifting toward 45° as 
soon as fi is off of the plateau. However, g(k) continues 
to peak at (ir, 0), and the product shifts off of (it, 0) more 
slowly as r is reduced. 

In conclusion by rather general symmetry arguments 
we have shown that d-wave superconductivity is a robust 
feature of the cuprates both hole and electron doped. We 
have also examined the evolution of the shape of the order 
parameter with doping and found a deviation of the or- 
der parameter angle dependence from simple cos 20-form 
similar to that measured in recent experiments (Fig. 4). 
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